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Abstract. We study the late-time evolution of the Universe where dark energy (DE) is 
presented by a barotropic fluid on top of cold dark matter (CDM). We also take into account 
the radiation content of the Universe. Here by the late stage of the evolution we refer to the 
epoch where CDM is already clustered into inhomogeneously distributed discrete structures 
(galaxies, groups and clusters of galaxies). Under this condition the mechanical approach is 
an adequate tool to study the Universe deep inside the cell of uniformity. More precisely, 
we study scalar perturbations of the FLRW metric due to inhonrogeneities of CDM as well 
as fluctuations of radiation and DE. For an arbitrary equation of state for DE we obtain a 
system of equations for the scalar perturbations within the mechanical approach. First, in 
the case of a constant DE equation of state parameter w, we demonstrate that our method 
singles out the cosmological constant as the only viable dark energy candidate. Then, we 
apply our approach to variable equation of state parameters in the form of three different 
linear parametrizations of w, e.g., the Chevallier-Polarski-Linder perfect fluid model. We 
conclude that all these models are incompatible with the theory of scalar perturbations in 
the late Universe. 
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1 Introduction 

Explaining the accelerated expansion of the late Universe is one of the greatest challenges 
in modern cosmology. The ACDM model is successfully described by the Planck data [1]. 
However, the nature of the cosmological constant is still unclear. In fact, there are about 120 
orders of magnitude between the observed value of the cosmological constant and the theo¬ 
retically expected value. The latter is related to the vacuum energy density. The cancellation 
mechanism that would reduce this large theoretical value of the cosmological constant is still 
a mystery [2], In addition, the ACDM model (as well as a lot of other dark energy models) 
faces the coincidence problem, i.e. why is the cosmological constant at present of the same 
order of magnitude as the dark matter energy density? To solve these problems, different dy¬ 
namical dark energy models were proposed. For example, such dark energy can be modeled 
by a barotropic perfect fluid (i.e. a fluid whose density is a function of its pressure only) with 
a proper equation of state (EoS). The Chevallier-Polarski-Linder (CPL) parametrization of 
the EoS [3, 4] is often used in literature (see, e.g., the recent papers [5-7]). In this approach 
the parameter of the EoS is a simple linear function with respect to the scale factor a of the 
Universe: w(a ) = wq ± wq(l — a/ao), where rcoy are free parameters of the model and a o 
is the present day value of the scale factor. The energy density £cpl of this fluid behaves 
as £cpl ~ a -3(i+«Jo+^i) a £ ear iy times when a <C a o, and when a is close to ao as well. 
Therefore, if rco,i are close to —1 and 0 respectively, then this barotropic fluid mimics the 
cosmological constant. Recently BOSS collaboration gave the constraints on the parameters 
of this EoS from the Planck+BAO+SN data as wq = —0.93 ±0.11 (la) and w\ = —0.2 ±0.4 
(la) [8]. In [9] an alternative form of the EoS was suggested following from the linear time 
parametrization w{t) = wq ± w\ (1 — t/to). Using the same data and method with [8], the 
authors found that the Planck±BAO±SN data predict for this model wq = —0.99±0.06 (la) 
and — 1 < wq ± w\ < —0.42 (2a). 

It is of interest to investigate these models from the point of their compatibility with 
other cosmological and astrophysical data, as well as with the theory of cosmological per¬ 
turbations, in particular, with the theory of scalar perturbations in the late Universe. At 
the late stage of the Universe evolution when inhomogeneities (such as galaxies and their 
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groups) are already formed, the hydrodynamic approach is inadequate. Here the mechanical 
approach [10, 11] is more appropriate. It works well inside the cell of uniformity [12] and 
provides us a good tool to investigate scalar perturbations for different cosmological mod¬ 
els (see, e.g., [13-15]). Therefore, it makes sense to study cosmological models filled with 
perfect fluids which have linear EoS parametrizations and investigate the compatibility of 
these models with the mechanical approach. Obviously, the above-mentioned forms of the 
EoS are among the simplest ones. In our paper, for completeness, we also consider a linear 
parametrization with respect to the redshift [16]: w(z) = wq — w\z. As a result, we show that 
all three considered models are incompatible with the theory of scalar perturbations in the 
late Universe. This means that such parametrizations can be used only for a limited period 
of time, starting, e.g., from last scattering till present. That is they follow from some more 
general EoS and for late times (in the future) they should be replaced by other expressions. 

The paper is structured as follows. In Sec. 2 we derive the basic equations for scalar 
perturbations in the framework of the mechanical approach in the case of a barotropic perfect 
fluid with an arbitrary EoS px = Px(^x)- Then, in Secs. 3, 4 and 5 we demonstrate that the 
enumerated DE EoS parametrizations linear in a, z and f, respectively, are incompatible with 
the scalar perturbations theory. The main results are briefly summarized in concluding Sec. 
6. In Appendix we consider a particular case of a barotropic perfect fluid with a constant 
EoS parameter. Here we investigate the conservation equations for scalar perturbations in 
the mechanical approach and demonstrate its self-consistency. 

2 Basic equations 
2.1 Background 

We consider the Universe at late stages of its evolution when galaxies and clusters of galaxies 
have already formed. At scales of the order of approximately 190 Mpc [12], being much 
larger than the characteristic distance between these inhomogeneities, the Universe is well 
described by the homogeneous and isotropic FLRW metric and the Friedmann equations 

-^2- — K (^O(CDM) + Dad + AxJ + A (2.1) 

and 

2 n' + n 2 + ic 

2 — -K (Prad + Px) + A > (2-2) 

a z 

where 7~L = a!/a = (da/dp)/a, k = 8irG]\r /c 4 (c is the speed of light and Gat is the Newto¬ 
nian gravitational constant) and K, = —1,0,+1 stands for open, flat and closed Universes, 
respectively. Conformal time p and synchronous time t are connected as cdt = adp. Further, 
T\ (cdm) is the average energy-momentum tensor of the pressureless (dust-like) matter. For 

such matter the energy density T^(cdm) = PcT 2 /® 3 is the only nonzero component, here 
p c = const is the average comoving rest mass density [10]. As usual, for radiation we have 
the EoS p rad = (l/3)e ra d- The Universe is also supposed to be filled with a barotropic perfect 
fluid with the EoS = Pxi^x) corresponding to DE. For generality, we also included the 
cosmological constant A into the above equations. Obviously, if the barotropic fluid is the 
reason of the late time acceleration of the Universe, then the cosmological constant is not 
necessary. It can be easily seen that the presence (or absence) of the cosmological constant 
does not affect the scalar perturbation analysis. For example, one can write down the Ein¬ 
stein equations with an upper index and a lower one. Then, the A term will correspond to 
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the contribution A 5 l k whose perturbations are zero at any order 1 . However, the dynamical 
behavior of the Universe (i.e. the form of the scale factor a) depends on A. 

From Eqs. (2.1) and (2.2) we can easily get the following auxiliary equation: 

^2 fa' — — K) = —K ( T 0 (CDM) + ^rad + + Prad + PXJ ■ (2-3) 


2.2 Scalar perturbations 


Deep inside the cell of uniformity the Universe is highly inhomogeneous. Here the mechanical 
approach is more adequate than the hydro dynamical one [10, 12]. In the framework of the 
mechanical approach galaxies, dwarf galaxies and clusters of galaxies (composed of baryonic 
and dark matter) are considered as separate compact objects. Moreover, at distances much 
greater than their characteristic sizes they can be described well as point-like matter sources. 
This is a generalization of the well-known astrophysical approach [18] (see §106) for the case of 
the dynamical cosmological background. Usually, the gravitational fields of galaxies are weak 
and their peculiar velocities are much smaller than the speed of light. These inhomogeneities 
together with fluctuations of other matter sources result in scalar perturbations of the FLRW 
metric [19, 20]. In the mechanical approach and in the conformal Newtonian gauge, the 
gravitational potential (i.e. the perturbation of the metric coefficient goo) rs_/ 0(1/c 2 ) 
satisfies the following system of equations (see [10, 12, 14] for details): 

A4» - 3 H(& + U®) + 3;/C$ = \^a 2 (dT 0 0 (CDM) + <fex + feradi + fc rad2 ) , (2.4) 

^(*' + W*) = o, (2.5) 

+ 3ft$' + (2 m' + H 2 )$ -K$ = i KCL 2 (dpx + dpradl + Sp rad2 ) , (2.6) 


where the Laplace operator A is defined with respect to the conformal spatial metric, and 
x@, (3 = 1,2,3, are the comoving spatial coordinates. From Eq. (2.5) we get immediately 
that 

®(U, r) = -^T, (2.7) 

Cra\Tj) 

where g?(r) is a function of all comoving spatial coordinates, and we have introduced c 2 in 
the denominator for convenience. In the surrounding of an inhomogeneity, the comoving 
potential <^(r) ~ 1/r [10, 12, 13], and the nonrelativistic gravitational potential ^(ry, r) ~ 
l/(ar) = 1/R, where R = ar is the physical distance. Hence, has the correct Newtonian 
limit near the inhomogeneities. 

Concerning the fluctuations of the matter sources, 4T ( /, CDM j is related to the fluctuation 
of the energy density of dust-like matter and has the form 2 [10] 


01 0 (CDM) 


5p c c 2 3 p c c 2 § 

Q ' Q 


( 2 . 8 ) 


x The conservation equation clearly shows that the A term can vary only in the case of interaction with 
other matter (see, e.g., [17]). However, this is out of the scope of our model. 

2 It can be easily seen from Eq. (2.23) in [10] that these fluctuations read (5To( CDM ) = 8p c c 2 /a 3 -f- 

(p c c 2 /a 3 )[34> + a$ 2 + /3® 3 + ...] where a and /3 are some constants. Taking into account Eq. (2.7), we 

get 4>/a 3 ~ 0(l/a 4 ), 4> 2 /a 3 ~ 0(l/a 5 ), 4> 3 /a 3 ~ 0(l/a 6 ) and so on. Since we consider linear perturbations, 

we should drop the terms 4> 2 ,4> 3 ,.... Therefore, in what follows, we will keep terms up to the order 0(l/a 4 ) 

inclusive. 
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where 5p c is the difference between the real and average comoving rest mass densities: 
5p c = Pc ~ ~Pc- We also split the fluctuations of radiation into two parts. The part la¬ 
beled by “radl” is caused by the inhomogeneities of dust-like matter (e.g., by galaxies and 
their groups), and the part labeled by “rad2” is related to fluctuations of the additional per¬ 
fect fluid (labeled by X). Of course, in the experiments we measure the total fluctuation of 
radiation. However, inhomogeneities associated with different matter sources can contribute 
to this total fluctuation separately. With respect to the additional perfect fluid, as we will see 
below, such contributions may take place directly via nonzero 5e ra ,d2 (cf. Eq. (2.22)). This 
happens if the inhomogeneities of the additional perfect fluid are distributed independently 
of the fluctuations of CDM [14]. However, there is also the possibility that the fluctuations of 
this perfect fluid nest on the inhomogeneities of CDM (such as galaxies, groups of galaxies, 
etc) screening the latter. In such models, <5e ra d2 = 0, <5ex ~ P and Eq. (2.17) (see below) is 
reduced (for physically consistent cases) to the Helmholtz equation [13]. Then, the fluctua¬ 
tions of the additional perfect fluid contribute to the fluctuations of radiation indirectly via 
the comoving gravitational potential <p(r) (cf. Eq. (2.10)) which is a solution of the corre¬ 
sponding Helmholtz equation. Given these two possible scenarios, we include 5e ra d2 hi °ur 
equations. Only a detailed analysis of a given model reveals which of these two scenarios is 
realized. 

For both of these contributions we have the same equations of state: 


^Pradl — (l/3)(5£radl j i(Pra<]2 — ( l/3)<l£rad2 ■ 


(2.9) 


We have shown in [12] that <5e r adi has the form 

<feradl=-^. (2.10) 

a 4 

It is worth noting that 5e r a d 2 should also behave as 1/a 4 (see Appendix). 

Let us now analyze Eqs. (2.4) and (2.6). Taking into account Eq. (2.7), we can rewrite 
them as follows: 


A4> + 3/C 4? — — Ata 2 + ^O-adl + <5£rad2^ , 

{H' - H 2 - K) 4> = )-Kd 2 ( 5px + 5p radl + 5p iad2 ) . 


( 2 . 11 ) 

( 2 . 12 ) 


Substituting Eqs. (2.8) and (2.10) into the right hand side (rhs) of Eq. (2.11), we obtain 


A4> + 3/C<f> = i kcl 2 


5p c c 2 

- n -b 0£x + 0£rad2 

a d 


(2.13) 


On the other hand, substituting (2.3) into the left hand side (lhs) of Eq. (2.12), we get the 
equation 

— ^^3-^ £ rad + e X + Prad + Px ^ ^ = $PX + ^Pradl + <^Prad2 • (2-14) 

Now, for consistency, we must take into account terms up to the order 0(l/a 4 ) inclusive. 
This is the accuracy of our investigations. Then, the terms E^dd? and p rad <l> of the order 
0(l/a 5 ) (cf. Eq. (2.7)) should be dropped: 

— ^^3-b £X + Px^j ^ = 3PX + ^Pradl + <5Prad2 ■ (2.15) 
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This equation can be rewritten in the form 


- (ex + Px) $ = dpx + gfe ra d 2 , 

where we have used Eqs. (2.9) and (2.10). 

Therefore, Eqs. (2.4) and (2.6) are reduced to the following system: 


(2.16) 


4 2 S 2 S 

a . KCra° c KC^a° _ . < 

Aip + 3 ICtp = -y-Spc H--— dsx H--—fe ra d2 , (2-17) 

— (ex + Px) ~J~ = 5px + o^ e rad 2 • (2-18) 

c z o 3 


Obviously, DE described by the barotropic perfect fluid must satisfy this system at any 
moment of the cosmic time in the late Universe. Such analysis should be performed for each 
EoS px = Px(ex) individually. 

For example, we can consider a perfect fluid with a constant parameter in the linear 

EoS: 

px = w£xi w / 0. (2-19) 

We exclude the case w = 0 since the CDM component was already taken into account in 
Eq. (2.1) 3 . From the conservation equation we have 


A 

£X ~ a 3 (l+w) ■ 

Further, for fluctuations of this matter we have 


( 2 . 20 ) 


5px = w5ex ■ 


( 2 . 21 ) 


Therefore, Eq. (2.18) gives 


p l + w _ 1 

0£X =-n-^X - — <)£rad2 

era w 6w 


( 2 . 22 ) 


Substituting this relation into (2.17), we arrive at the following equation: 


A tp + 3 Kp 



nip 1 + w 
2 w 


a 2 ex + 


kc 2 o 3 
2 



(2.23) 


The lhs of this equation does not depend on time 4 . Therefore, the same should hold true 
for its rhs. Because <5e ra d2 ~ l/ fl4 > w e must demand that <fe ra d 2 = 0 (there is no possibility 
that the term with Ex and the term with <5e ra d2 cancel each other). Therefore, there are only 
two values of w which do not contradict this equation. They are w = —1 and w = —1/3. 
The former case reduces to the standard ACDM model considered in [10]. The latter case 
w = —1/3 was considered in detail in [13] in the absence of radiation. The same conclusions 

3 In addition, if w = 0: (i) we get from Eq. (2.18) that fc ra d 2 = —£x[3</?/(c 2 a)] = —3 A<p/(c 2 a 4 ) which is 
similar to Eq. (2.10), and (ii) Eq. (2.17) implies &x = B/a 3 — & ra d 2 in full analogy with Eq. (2.8) (for 
B = 5p c xc 2 ). Obviously, we can combine both of these dust components into a single one. Therefore, we 
would be back to the model with CDM and a A term (in the background equations (2.1) and (2.2)). 

4 We would like to remind that the quantities ip and 5p c are the comoving ones [10]. Therefore, within the 
adopted accuracy when both nonrelativistic and weak field limits are applied, they do not depend explicitly 
on time [12]. 
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with respect to the model with constant w can be obtained from the conservation equations 
(see Appendix). Obviously, the dark energy component is smooth and homogeneous in the 
background equations where matter/energy is considered in the averaged form. However, 
it does not have to be smooth and homogeneous a priori. The conclusion of the presence 
or absence of dark energy fluctuations follows from the analysis of Einstein equations for 
perturbations. For example, for the ACDM model (i.e. w = — 1) fluctuations are absent, as 
it is easily seen from Eq. (2.22) where <5e ra d2 = 0 (see also [10]), so dark energy represented 
by the cosmological constant A is truly homogeneous, while for the model with w = —1/3, 
fluctuations are nonzero [13], so there are inhomogeneities in the dark energy distribution 
over the Universe. It is clear that in the case of the inhomogeneous distribution, the average 
value of fluctuations should be equal to zero: dex = 0 (see, e.g., [13] where it is shown 
explicitly). 

We next investigate the models with dynamical w. 


3 CPL model 

First, we consider the CPL model (i.e. X = CPL) and the corresponding EoS [3, 4]: 


Pcpl = w(a)e cpl , w(a) = w 0 + w 1 [l- 


« o 


(3.1) 


where ao represents the scale factor of the Universe at the present time. We consider this 
EoS to be valid for a definite period of time in the past (e.g., from last scattering [4]). We 
also suppose that this EoS is still valid for some period of time in the future. 

From the conservation equation in the form 


d(s cpl 0 - 3 ) +PcPL^( a3 ) — 0 
we easily get the well known formula 

e CPL = ^ a -3(l+^0+^l) e 3^ia/a 0 ^ 

where A is a constant of integration. Therefore, 


Pcpl + £cpl = Aa~ 3 ( 1+Wo+Wl ^ e 3wia - 


/a 0 


l + t«0 + Wl 1- 

a 0 


(3.2) 

(3.3) 

(3.4) 


Additionally, with the help of Eq. (3.3), we obtain the expression for the fluctuation of the 
CPL fluid pressure: 


fipcPL = 


d (wscpl) 

decPL 


&TCPL = ( w - 


fecPL 


W 0 + W! 1-- - 


ao 


WjgCPL da \ 

«o de cpl / 

1 W\ 

3»i - ^ (! + w 0 + wq)_ 


decPL • 


(3.5) 


Therefore, Eq. (2.18) reads 


A( P_- 3 


_ _ a -3(l+uio+uii)-lg3uiia/a 0 

C 2 


1 + Wo + W\ 1 — 


wq T w\ [ 1 —— I — 


a o 


w 1 


3ioi - ^ (l + w 0 + wi)_ 


a o 


deCPL + 2<ferad2 


(3.6) 
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and for fecPL we have 


a^ dt = A v^sa+wn+w^-l ^B (1 + Wo + w± (1 — /3)) 2 

c 2 w 0 + w 1 (l-2p/3) + (w 0 + w l (l-P)) 2 

_ 1^ £ _ 1 + w 0 + wi (1 - /3) _ 

3 rad2 w 0 + wi(l- 2/3/3) + (w 0 + u>i (1 - /3)) 2 ’ 


(3.7) 


where for convenience we have introduced the notation /3 = a/a q. Consequently, Eq. (2.17) 
takes the form 


Al p + 3 K V = ^-S Pc - ^ a -3(l + ™ + ».) + 2 e 3„^_(1 + WQ + W! (!-<?)) 


w 0 + Wi (1 - 2/3/3) + (tco + Wi (1 - P)Y 


+ 


nc 2 a 3 


<fe r ad2 


1 

1- 


1 + w 0 + w\ (1 - /3) 


3 w 0 + w\ (1 - 2/3/3) + (w 0 + w\ (1 - /3)) 


(3.8) 


The particular case w± = 0, i.e. w = wq = const, was considered at the very end of 
the previous section. Now, we turn to the general case w i / 0 and consider Eq. (3.8). Here 
the lhs and the term containing 5p c are independent of time. Therefore, either each of two 
remaining expressions are also independent of time (within our accuracy 0(l/a), bearing in 
mind that the rhs of (3.8) has been multiplied by a 3 ), or they must cancel each other at any 
arbitrary moment of time. Can this be the case? The second assumption that the second and 
third terms on the rhs of Eq. (3.8) can cancel each other at any arbitrary moment of time 
does not work because of the presence of the exponential (with respect to the scale factor a) 
function exp(3uq/3) which is linearly independent from the power functions. It can be easily 
seen that the first assumption does not work as well in the case w± Y 1 0. For example, we 
can eliminate the last term in Eq. (3.8) by putting <3e ra d2 = 0. However, the second term on 
the rhs of this equation still depends on time. The simple analysis shows that, in contrast 
to Eq. (2.23), we cannot select the parameters rco,i of the model in such a way that this 
term becomes independent of time. Already here, we can see a potential problem of the 
model. However, this term can still be of the order o(l/a) and then we can drop it within the 
accuracy of our approach. The problem is that we cannot expand the exponential prefactor 
into series with respect to the scale factor which can be arbitrary large. Therefore, we need 
to estimate the second term taking into account the exponential function. Obviously, this 
term is exponentially small for w\ < 0 and (3 l 5 . However, we cannot drop the second 

term simply on the basis of this argument. To show it, we consider the following ratio: 


—3(l+wo+wi)+2 3wi/3 _ 

Li O 9 C 

2 AvC^Q" 


1 


;2/^„-3(1+u>o+uji)+2 „3iui/3 

2 


^radl 

1 


ECPL 


,3£c£ 3 p c c 2 /a 3 


(3.9) 


Now, we assume that CPL matter is responsible for the late time acceleration of the Universe. 
According to the recent observations, it contributes approximately 70 % to the total energy 
density of the Universe, while CDM + baryons give approximately the remaining 30%. 
Hence, at the present moment the ratio (3.9) is approximately equal to 0.78. If we consider 
for wo and w\ the values —0.93 and —0.2 [8], respectively, then at the present moment (i.e. 

5 For positive values of wi > 0, the second term is exponentially divergent when a —> oo. Hence, there is 
no need to consider such values since our model must be valid for arbitrary large values of the scale factor a. 
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for (5 ~ 1) we get that the second term amounts 3% from the radiation which we take 
into account in our analysis. This amount is not negligibly small and we cannot neglect 
such a time dependent term. Therefore, the considered CPL barotropic fluid with rai / 0 
contradicts the theory of scalar perturbations in the late Universe. 

4 Dark energy with an EoS parameter linear in z 

Now, we consider a model where the parameter of the EoS for DE is a linear function of the 
redshift z [16]: 

p z = w(z)e z , w(z) = w 0 — w\z = w 0 + w± ^1 - —^ . (4.1) 

For this perfect fluid we get 

= Aa~ 3{1+W0+Wl) e ~ 3wiao/a . (4.2) 

By performing an analysis similar to the one carried on the previous section, we conclude 
that this model also contradicts the theory of scalar perturbations. Obviously, it happens 
because of the exponential factor in Eq. (4.2) which leads to an equation for the comoving 
gravitational potential similar to (3.8) where the last two terms on the rhs do not cancel each 
other. 


5 Dark energy with an EoS parameter linear in t 


In the two previous parametrizations of the EoS w is a divergent function either for a —>■ +oo 
or for a —> 0. In the recent work [9] a different parametrization for w which is based on a 
linear dependence on time t was proposed: 


w(t) = Wo + W\ 1 - 

to 


leading to the following expression in terms of the scale factor a: 

2(1 + w 0 + wi) 


w t (a) = w 0 + wi 


1 - 


wi + (2 + 2u) 0 + w 1 )/3- 3 ( 1+wo+wi y 2 _ 
—wi (w 0 + w\ + 2 ) + (2 w 0 + w\ + 2)(rt> 0 + u>i) 

wi(3 3uj / 2 + 2wq + w\ + 2 


(5.1) 


(5.2) 


It was shown in [9] that if DE is parametrized by the EoS (5.2) then: (i) data from Planck, 
BAO and SN constrain much better the free parameters of this model (i.e. u>o and w±) as 
compared with the CPL parametrization, and (ii) if (tco, w\) = (—1, 1/3) in Eq. (5.2), i.e. 
the DE EoS parameter varies between —2/3 in the early Universe and —1 at present, it fits the 
data slightly better than a pure A term. The relation (5.2) was obtained under the assumption 
w ^4 0, where w = wq + w\ + 1. More precisely, 0 < w\ < 2w. In general, we can analytically 
continue this expression to negative values of w, keeping the condition w ^ 0. As we will see 
below, physical quantities (e.g., the energy density) have the same asymptotic behavior at 
big values of the scale factor a for both positive and negative w. The EoS parameter (5.2) 
goes to a constant value in the limit a — > +oo: wt(a) —> —(w + l),w > 0 and wt(a) —> 
— (|u)| + 1), w < 0. In the case w\ = 0 we naturally have wt(a) = wq = const which was 
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considered at the end of Sec. 2. We note that the direct substitution of w = 0 into Eq. (5.2) 
results in an uncertainty of the type 0/0. One can check that wt(a) —> as 

wq —> — 1 — w \. This is not a rational fraction and we again will have a problem (except for 
the trivial case wq = —1) similar to the one described in two previous sections. Therefore, 
we disregard this case. 

The conservation equation for the considered fluid reads 

d{E t p 3 ) + w t £ t d(P 3 ) = 0 . (5.3) 


Now, with the EoS parameter (5.2), we get the energy density 


— £t 0 


f wifi 3 ™! 2 - w i 

1 + m 


4 




(5.4) 


where £t 0 is an integration constant. It can be easily seen that for both positive and negative 
w the energy density behaves asymptotically as Et p3\w\ -(-qq for j3 —+oo. In the case 

w\ = 0 this formula is reduced to the familiar expression: Et = Et 0 /3~ 3 ( 1+W °\ 

Let us turn now to the scalar perturbations for this model. Taking into account Eq. (5.4), 
we get 
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With the help of Eq. (2.18) we obtain the expression for b£t- 
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Then, Eq. (2.17) takes the form 
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It can be easily seen that this equation reduces to Eq. (2.23) for the particular value w\ = 0 
and, consequently, w = wq. Let us now consider the general case w\ / 0. According to 


- 9 - 































the accuracy of our approach, we should take into account all terms up to the order 0{l/a) 
inclusive (bearing in mind that the rhs of (5.8) has been multiplied by a 3 ). It should take 
place for all times in the late Universe including the far future (3 1. In this limit it is very 

convenient to expand the expressions under consideration in powers of the scale factor. We 
see that, because the first line in Eq. (5.8) is independent of time (or of the scale factor), all 
terms of the order 0(l/a n ), n < 1 in the next lines must be absent or they must cancel each 
other. The term in the second line behaves asymptotically as rsj /3 3 H+2, w h en /3 1, i.e. its 

asymptotic behavior does not depend on the sign of w. The last term behaves asymptotically 
as ~ /3 3 de ra d2, if /3 3> 1, again for both positive and negative w. Therefore, this term has the 
asymptotic form ~ 1/(3 in the presence of <5e ra d2 ~ 1/a 4 or it equals identically to zero for 
<ferad 2 = 0. Hence, first, in the presence of fe r ad 2 the second line and the last term cannot 
cancel each other, and, second, in the absence of such radiation the time dependent second 
line cannot be dropped from the equation. This result is in contradiction. We conclude 
that the problem has arisen because of the power-law growth of the energy density in the 
expanding Universe: Et ~ /3 3 ^ when /3 3> 1. Therefore, the model considered in this section 
contradicts as well the theory of scalar perturbations in the late Universe. 

6 Conclusions 

In our paper we have studied the late Universe filled with dust-like matter and a barotropic 
perfect fluid which can play the role of dark energy. We have considered the radiation content 
of the Universe as well. More precisely, we have concentrated on scales much smaller than the 
cell of uniformity size which is approximately 190 Mpc [12]. At such distances our Universe 
is highly inhomogeneous and the averaged Friedmann approach does not work here because 
we need to take into account inhomogeneities in the form of galaxies, groups and clusters 
of galaxies. It is also a natural assumption that radiation and the barotropic perfect fluid 
fluctuate around the average values. Therefore, these fluctuations as well as the dust-like 
matter inhomogeneities perturb the FLRW metric. To analyze these perturbations inside 
the cell of uniformity, we need to use the mechanical approach [10-12]. An important feature 
of this approach is that it provides an optimal tool to study the compatibility of different 
cosmological models with observations (see, e.g., [13, 14]). 

In the present paper we have firstly derived the equations for scalar perturbations within 
the mechanical approach in the case of a barotropic fluid with an arbitrary EoS px = Px(ex) 
(on top of the dust-like matter and radiation content of the Universe). Our master equations 
are Eqs. (2.17) and (2.18). First, we have shown how these equations work in the simplest 
case of a constant EoS parameter w for dark energy. Our analysis demonstrates that the only 
two values of w, namely w = —1 and w = —1/3, are acceptable. Therefore, in the case of 
constant w, the cosmological constant (i.e. w = —1) is the only viable dark energy candidate 
(since a perfect fluid with w = —1/3 does not accelerate the Universe). 

Then we have applied these equations to variable EoS parameters. Since it is impossible 
to handle this problem in full generality, we consider three linear parametrizations of w: CPL 
(i.e., the parameter of the EoS is linear in the scale factor a), linear in the redshift z and linear 
in time t. Our analysis has shown that all these three models contradict the theory of scalar 
perturbations in the late Universe. This occurs despite the fact that these models can be in 
sufficient agreement with the observed data obtained from Planck, BAO and SN [1, 9, 16]. 
This means that these EoS are not the fundamental ones and they can be used only for a 
limited period of time, starting, e.g., from the time of last scattering and up to the present 
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moment. During this period they approximate/mimic more or less successfully the unknown, 
’’real” EoS for DE. In the future Universe they should be replaced by other expressions. 
Our investigations show that the most simple and usual constituents of the Universe such 
as the cosmological constant, pressureless dark and baryonic matter and radiation are in 
full agreement with the theory of scalar perturbations. Interestingly, a perfect fluid with a 
constant EoS parameter w = —1/3 fits this scheme also very well [13]. 

We highlight that our conclusion on the incompatibility of perfect fluids which have 
variable EoS parameters with the theory of scalar perturbations in the late Universe applies 
only to the considered cases of linear parametrizations of w. In our study, we assumed that 
these parametrizations take place up to very big values of the scale factor a. This is a rather 
strong restriction. Therefore, the cases of variable w different from the linear parametrizations 
require a separate investigation. For example, it was shown in [14] that exotic quark-gluon 
nuggets suitable to describe the dark sides of the Universe [21, 22], if they exist, can also be 
in agreement with the theory of scalar perturbations while being characterized by variable 
w. 
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Appendix: Adiabatic perturbations of perfect fluids with constant EoS pa¬ 
rameters 

In this appendix we obtain the fluctuations of the energy density for perfect fluids with a 
constant EoS parameter as a function of the scale factor a. For this purpose we consider the 
conservation equation (see, e.g., Eq. (2.74) in [20]) 

de'x + 3— (Sex + Spx) + (ex + Px) (Av — 3H/) = 0, (Al) 

a 

where A is the Laplace operator with respect to the comoving coordinates and v is the 
comoving peculiar velocity potential. We remark that we consider the fluctuations of perfect 
fluids arising due to inhomogeneities of dust-like matter (e.g., galaxies). These fluctuations 
also form their own inhomogeneities. Then the velocity potential v characterizes dynamics of 
displacement of such inhonrogeneities. For example, in the case of radiation, v corresponds 
not to the speed of light but to the speed of displacement of photon concentrations. In the late 
Universe these displacement speeds are nonrelativistic and give much less contributions than 
$ [10, 12], Therefore, we drop the term containing v in Eq. (Al). The gravitational potential 
is given by the formula (2.7): = p/(c 2 a ). We consider perfect fluids with a constant EoS 

parameter w : px = wsx- Hence, the conservation equation gives ex = A/a^ l+w \ For such 
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fluids their adiabatic fluctuations also satisfy the same EoS: 5px = w5ex■ Then the general 
solution of Eq. (Al) may be written in the form 

c- A 3y(l +w)A 

x 0 3(i+«,) ^ c 2 a 3(i+«,)+i ’ 

where A is a constant of integration. This expression shows that in the case of the negative 
parameter w < 0 both terms in the rhs satisfy the adopted accuracy because 3(l + u;) + l < 4. 
However, the Einstein equations (in our case, in the form of Eqs. (2.17) and (2.18)) may 
impose additional restrictions on the form of 5ex- For example, the agreement with Eq. 
(2.22) where <5e r ad2 = 0 results in additional conditions: A = 0, —(1 + w)/w = 3(1 + w). The 
latter condition holds true only for w = — 1 and w = —1/3. Therefore, 

fe -i/3 = ' ( A3 ) 

The conclusion that the value w = —1/3 is a particular one was also made from different 
reasoning in papers by Fulvio Melia (see, e.g., [23] and the references therein). It can be also 
easily seen that in the case w = 0 the formula (A2) exactly reproduces Eq. (2.8) if we put 
A = 5p c c 2 , A = p c c 2 . In the case of radiation w = 1/3 the second term in (A2) is of the 
order 0(l/a 5 ) and should be dropped. Then, we get 

fcrad ~ 1/a 4 • (A4) 

It is worth mentioning that in the general case of interacting fluids Eq. (Al) must 
contain their total energy densities and pressures. This case is beyond the scope of our paper 
where we concentrated exclusively on the standard perfect fluids without energy exchange 
between them. 

There is also one more conservation equation (see, e.g., (2.75) in [20]). In our approach 
(that is when neglecting the peculiar velocities) this equation reads 

+ + Px*) $ ] = 0 • ( A5 ) 

i 

Let us consider a model with four perfect fluids: dark energy represented by the cosmological 
constant, CDM, radiation and the w = —1/3 fluid. Then, (A5) takes the following form: 

5pA + (£a + Pa) + ^pcdm + (ecdm + pcdm) 

+ ^Prad + (e r ad + Prad) ^ + $P-l/3 + (^—1/3 + P- 1/3) = 0 . (A6) 

Because 6p\ = 0 and E\ = — p\, the contribution from the A-term disappears. Now, taking 
into account that dpcDM = Pcdm = 0, 5p rad = (l/3)<5e rad and (e rad +p rad ) 4* ~ 0(l/a 5 ), we 
get 

ferad = -3ffCDM^ =-, (A7) 

where £cdm = p c c 2 /a 3 . This expression exactly coincides with the formula (2.10) in Sec. 2. 
Finally, keeping in mind that p_ 1 / 3 = —(l/3)e_ 1 / 3 and <5p_i/ 3 = — (l/3)5e_ 1 / 3 , we arrive at 
Eq. (A3). Therefore, we have shown the self-consistency of our approach. 
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